A new many-body wave function for BCS-BEC crossover in Fermi gases 
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We present a new many body formalism for BCS-BEC crossover, which represents a modification of the 
BCS-Leggett ground state to include 4-fermion, and higher correlations. In the BEC regime, we show how our 
approach contains the Petrov et al 4-fermion behavior and associated scattering length add at short distances, 
and secondly reduces to composite-boson Bogoliubov physics at long distances. It reproduces the Lee- Yang 
term, whose numerical value is also fixed by add- We have also examined the next term beyond the Lee- 
Yang correction in a phenomenological fashion, building on cloud size data and collective mode experiments. 
However, one has to view this phenomenological analysis with some caution since experiments are in a state of 
flux and are performed close to unitarity. 
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There have been exciting new developments in the field 
of trapped fermionic atoms 01 01 in which applied mag- 
netic fields are used to effect a transition from BCS superfluid- 
ity to Bose Einstein condensation (BEC). The simplest mean 
field wave function 0, 01 which describes these phenomena 
at T = is, effectively, the BCS wave function |"0bcs) with 
self consistently determined (fermionic) chemical potential in 
the presence of pairing interactions of arbitrary strength. This 
has met with some experimental success |6] but it is clearly 
lacking in the sense that it does not incorporate the exact two- 
molecule (4-fermion) scattering physics |7]. In this mean field 
theory, the 4-fermion interaction, parametrized by the dimer- 
dimer scattering length ctdd> is overestimated by a factor of 3, 
leading to a similar overestimate in the cloud size |8]. Con- 
densate depletion is also absent in this lowest order mean-field 
theory. 

This earlier work raises a fundamental question which we 
address here. How does one arrive at a many body wave func- 
tion which incorporates the few body constraints and what are 
its implications for experiment. While, there have been efforts 
[loll to go beyond the leading order mean field theory in 
the literature, there has never been a demonstrated consistency 
between the exact few body constraints and N-body behavior. 
Our work is the first to make these two points of view compat- 
ible and in the process, introduces a new perturbation scheme 
which replaces the usual Feynmann diagrammatic expansion. 
Here we derive an equation of state for molecular bosons at 
the level of the Bogoliubov approximation, and, thereby set 
up a basis from which one can compute various properties of 
the gas. Of particular interest are those which are associated 
with the underlying fermionic physics, which are absent in a 
conventional BEC of true bosons. 

The parameters appearing in our wave function are de- 
termined directly by an energy minimization procedure; we 
stress that there are no fitting parameters. The equation of 
state reduces to the well known Lee- Yang like formula flUl . 
but with the scattering length add determined by the exact 4- 
fermion physics 01 . Finally, our wave function appears com- 
patible with the general features of experimental data on two 



different experiments: cloud sizes |8] and collective modes 

O. 

We begin by summarizing our results for the many body 
wave function and its consequences before providing the de- 
tailed methodology we employ. Additional details are given 
in the unpublished supporting material (Ref. 1 13]) which we 
have made available. One may rewrite the mean field wave 
function |4j,l3l as 



bcs) oc exp 



a k c kT c -k| 



|o>, 



(1) 



where destroys fermions with momentum Kk and spin a 
0, and |0) is the particle vacuum. Written in this way, the 
BCS wave function is seen to be the "composite-boson" ver- 
sion of the Gross-Pitaevskii (GP) wave function of conven- 
tional BEC. 

It seems natural to presume that corrections to this GP-like 
wave function should capture some of the physics of the con- 
ventional Bogoliubov approach which provides an improve- 
ment over the GP wave function, appropriate to bosonic sys- 
tems. We may rewrite the corresponding Bogoliubov wave 
function as 



Bogoliubov/ 



exp 



q>0 



10} , (2) 



where 6^ is the creation operator of bosons with momentum 
hq which, like b\, must be replaced by a bilinear form of 
fermion creation operators. In this way, one sees that four 
fermion terms are required to go beyond the simple BCS-like 
mean-field theory. 

On this basis we return to Eq. Q and write the natural cor- 
rection terms. As will be shown below, to capture the physics 
of the deep BEC region, 6-fermion as well as 4-fermion cor- 
rections are needed. This goes beyond earlier work |9]. More 
generally one can contemplate a many body wave function of 
the form 
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where each K; represents a shorthand notation for k^Ci, and 
— K refers to the fact that both the momentum and the spin are 
reversed. There are constraints imposed on the coefficients a, 
(3, 7, etc, which can be readily formulated 11311 . although we 
will not do so here. 

We show below how to compute physical properties in 
terms of these wave function parameters. We do this in a fash- 
ion which is also consistent with Ref. 0]. Thus, to determine, 
for example, add, the wave function is adjusted to minimize 
the energy, at a given fixed number of fermions. For general 
values of kp a this procedure is very difficult to implement. 
However, in the low-density regime, given more precisely by 
n d a dd ^ 1' higher-order correlations in the wave function 
are much weaker, and we are able to construct a systematic 
perturbation theory. For this reason, our theory is mostly lim- 
ited to the deep BEC regime of the crossover (rida dd <C 1). 
This is different from some other approaches which produce 
various smooth curves of the ground state energy in the whole 
BCS-BEC crossover 16]. Despite this limitation, phys- 

ical observables (such as the equation of state) take the form 
of well-controlled low-density expansions, and the functional 
forms, as well as the coefficients of such expansions can be 
exactly determined within our theory, whereas the alternative 
theories llfi HTslflfill lack such exactness in the BEC limit. 

The ground state expectation value of any physical observ- 
able O is 

(0|TOexp(qt +/? t + ... + a + /? + ... ) |o ) 
[ ' (0|T exp(at + (}t + . + a + + . .. )|0) ' 1 ; 

where a, /3, etc, are the terms in the exponent of the right 
side of Eq. (0, and T is a "time-ordering" operator defined 
by assigning the fermion creation operators in a, 0, etc the 
earliest "time", the annihilation operators in a', 0\ etc the 
latest "time". Finally, those times appearing in O enter as 
intermediate "times" increasing monotonically from right to 
left. A standard Feynman diagrammatic method is then used 
to evaluate Eq. @. It should be noted that, in our present ap- 
proach, the terms in the wave function (instead of the interac- 
tion) serve as vertices, in contrast to the textbook applications 
of Feynman diagrams. Our approach is suited to the bosonic 
regime, where the attractive interaction between fermions is 
so strong that they form bound states. 

One can contemplate very general interactions between two 
fermions which lead to these bound states. Under these cir- 
cumstances, it follows from Eq. that the ground state en- 



ergy is given by 

E/Q, = n d E d + 2Trh 2 n d add/m d 

x {l + [128/(15^)] + ■••}, (5) 

which is similar to that of point-like bosons full . Here f2 
is the volume, rid — n/2 is the number density of dimers, 
nid = 2m is the mass of the dimer, m is the mass of each 
fermion, E d °^ the energy of a zero-momentum dimer in vac- 
uum, and add is the self consistently determined dimer-dimer 
scattering length. It should be stressed that this equation is 
very general; it applies even when the fermionic interaction 
is not a contact interaction. When the interaction has a finite 
range (comparable to the fermionic scattering length a), the 
famous calculation of Petrov et al in Ref. |7] needs to be cor- 
rected and add/ a 7^ 0.6. It is the actual add, determined from 
the solution to the quantum 4-fermion problem for the low- 
energy scattering of two molecules [see Eq. l|9} below], that 
must be used in Eq. (|3J- When the interaction is of very short 
range (<§; a) then the add in our theory reduces to that in Ref. 
10. 

An equation of state associated with the first two terms in 
Eq. (|5} (between the curly brackets) was first proposed in Ref. 
11711 . based on the assumption that fermionic dimers were gov- 
erned by the same equation of state as point-like bosons. It 
is of considerable interest that we have confirmed this equa- 
tion of state in the context of a fermionic system. Also note 
that this equation of state [Eq. (0] contains ndd dd < 1 as 
a small parameter flill . While one might be inclined to con- 
sider a power series in kpa, in this low density regime, where 
kp cc n 1 / 3 is the formal Fermi wave vector, we have shown 
that this is incorrect. When all fermions form tightly bound 
pairs, there is no Fermi sphere, and no natural expansion in 
terms of n 1 / 3 exists. 

We now provide some details for our calculation by sketch- 
ing how this wave function is determined in the BEC regime. 
If ndadd ^ 1' higher correlation terms in the exponent in 
Eq. (0 are progressively less important. Their orders of mag- 
nitude are derived, determining which coefficients we should 
retain, at any given order in the expansion of the equation of 
state. At the leading order, only «k is important, and the two- 
body Schrodinger equation is reproduced upon minimization 
of E — fiN (N is the number of fermions). This wave func- 
tion corresponds to a collection of dimers. At the next order, 
enters and is adjusted to minimize energy. It is consequently 
found that satisfies the 4-body Schrodinger equation for the 
low-energy scattering of two molecules at the leading order. 
We thus use the solution to the latter equation to construct 
approximately [the result is shown in Eq. (|8) below]. In the 
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process ok is computed to a higher order, to minimize energy 
at this higher level of approximation. Summing up all the con- 
tributions to energy up to this level (using the results for a and 
j3 that one just obtained), one deduces a mean-field interaction 
between dimers; note that here one derives the correct add in 
accordance with Ref. 1 7] in such interaction term. 

At the third order, 7K 1 +q/4,-,K 4 +q/4,k'-q/2T ,-k'-q/2| 
now appears. This order corresponds to the usual Lee- Yang 
term 1 1 1 ] for point bosons, although here the derivation is for 
the molecular boson case. Similarly, ok and j3 are recom- 
puted to yet a higher order of accuracy. Corrections to the 
latter are now predominantly at small q [~ O^ndadd)] ■ 

At this third order, the function (3 in Eq. unifies the Bo- 
goliubov theory of composite bosons and the exact 4-fermion 
calculations of Ref. [7]. We show this as follows by deduc- 
ing an asymptotic formula from Bogoliubov theory for com- 
posite bosons. We introduce o q = <^kCq/ 2+k |C q , 2 _ k ^, 
where 0k is the internal wave function of the dimer, satisfy- 
ing X)kl^k| 2 = 1. With this substitution, it follows that 

/3q/2+kT,q/2-kJ.,-q/2+k'T,-q/2-k'l ~ ^q^k^k'- (6) 

In fact, by strictly applying the energy minimization princi- 
ple alone, without resorting to the picture of bosons, we are 
able to derive the same formula as above (under the condition 
q -C 1/Vd and q <C |k — k'|, where r<j ~ a is the molecular 
size), and our formula for x q coincides, at the leading order, 
with that of a traditional Bose gas full whose scattering length 
equals a d d- For l/r d > q > ^n d a d d and |k - k'| > q, our 
theory yields the standard asymptotic formula for x q to obtain 

-iim d add, , 

Pq/2+kT,q/2-kl,-q/2+k'T,-q/2-k'l — 2 9^9^' 

(7) 

plus higher order corrections. 

Now let us compute the related behavior which derives 
from the nonperturbative four-body physics of Ref. |7] . When 
q and |k — k'| are both ^> ^JridCidd, it can be shown that 

/5 q /2+kT,q/2-k|,-q/2+k'T,-q/2-k'J, 

= A^ q /2+kT,q/2-kl,-q/2+k'T,-q/2-k'i 

plus higher order corrections. Here Nd = N/2, and <p^ is the 
momentum representation of the four-fermion wave function, 
describing the low-energy scattering of two dimers. 

In coordinate space, such a four-fermion wave function 
takes the asymptotic form oc (1 — add/Rdd)<f>( r )<t>( r '), wnen 
the distance between dimers Rdd ^> I'd- Translating this long- 
distance boundary condition into momentum space, we get: 

^q/2+kT,q/2-kJ.,-q/2+k'T,-q/2-k'J. 

« (* q ,o - ^r)^*, i « Uu. (9) 

where |k — k'| ^> q, and ft is the volume. 



Note that Eq. is valid in the phase space region in which 
q <C l/rd, |k — k'|, while Eq. (jSJl is valid in the region 
q, |k — k'| ^> y/ridddd- Now consider the region of over- 
lap, in which l/rd ^ 1 \fndOd2 an d |k — k'| 3> q. In this 
region, both the composite-boson formula [Eq. ©J and the 4- 
fermion formula [Eq. 0] yield exactly the same result, Eq. Q 
above. Thus the conventional picture of many bosons (whose 
condensate suffers from a quantum depletion, as shown by the 
function a; q ) and the nonperturbative four-fermion effects are 
merged into a single picture. 

Finally, we turn to experimental implications of this ap- 
proach closer to unitarity. Here we introduce higher order 
terms in the equation of state [Eq. in a phenomenolog- 
ical manner. Our ansatz, for which we have some micro- 
scopic support , builds on earlier exact results fl9l Eoil for 
non-composite, or true bosons. Thus we presume that 

E/fl = n d E d ^ + 2irh 2 n d add/md 

x (1 + 4.8144V? + 19. 6539 ff In g + Cg) (10) 

plus higher order terms which we shall omit below. Here g = 
n d a dd i s m e B° se g as parameter. The coefficient in front of 
g In g was first found by Wu for the hard sphere Bose gas fiUl . 
and was later shown to be valid for all gases of identical non- 
composite bosons. 

The phenomenological term containing the unknown pa- 
rameter C is roughly of the same order of magnitude as the 
Wu term. Physically, the C-term must involve a three-dimer 
scattering process. This corresponds to a quantum six-body 
problem, making C very difficult to determine directly from 
theory. For a dilute gas of bosonic atoms with large posi- 
tive scattering length, a quantum three-body calculation l2lll 
yields C w 141. However, this latter result is not applicable 
to the BEC-BCS crossover system we study here. 

To proceed, we calculate C from currently available experi- 
mental data, on the radial breathing mode frequencies 1 12[ and 
axial cloud sizes |8] of cigar-shaped two-component Fermi 
gases of 6 Li atoms in harmonic traps. The atomic scattering 
length a is taken from Ref. 1 22] where it is found to diverge at 
B=834G. We presume a dd = 0.6a H . We start from Eq. JTOb . 
and evaluate the axial size (using the local density approxima- 
tion for the density profile) and the radial breathing mode fre- 
quency (using the scaling ansatz of Ref. 1 23 ] ). The parameter 
C is then adjusted to best reflect the general features of the 
data. In this way we arrive at a phenomenological equation 
of state which extrapolates from deep in the BEC (where it is 
well controlled) into the range of kpa accessible by experi- 
ment. 

In Fig.Q we plot the normalized cloud size £ as a function 
of B. In Fig. [2] we plot the radial breathing mode frequency 
as a function of 1 /kpa. In both figures, kp is the Fermi wave 
vector of the noninteracting Fermi gas at the trap center. A 
value C = 40 reproduces the trends of the first generation 
radial mode frequencies IT2I1 and is also not inconsistent with 
the cloud size data |8]. From the cloud size data, we observe 
that there may be some uncertainty in C of the order ±10, 
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FIG. 1: (color online). Axial cloud size <T ax iai normalized by that of 
the noninteracting Fermi gas (with the same number of atoms), of, 
for a cigar-shaped trapped Fermi gas of 6 Li atoms. Error bars: data in 
El. Solid line: Eq. (TOJ with C = 40. Dashed line: C = 50. Dotted 
line: C — 30. Sparse dotted line: C = 141, as in an atomic BEC 
with large scattering length |21], Dot-dashed line: bosonic mean- 
field approximation, as in 1 8] . In the inset are plotted the differences 
between the various cloud size curves and the bosonic mean-field 
curve. 
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FIG. 2: (color online). Radial breathing mode frequency (normal- 
ized by the radial trap frequency) versus 1/kFa. Error bars: data 
from 03]. Solid line: Eq. (l0| with C = 40. Dashed line: C = 50. 
Dotted line: C = 30. Sparse dotted line: C = 141, as in an atomic 
BEC with large scattering length [21]. Dot-dashed line: bosonic 
mean-field approximation. Dot-dot-dashed line: Monte-Carlo sim- 
ulations 1 24]. The solid line is extended to indicate the likely con- 
tinuous curve consistent with data, however there is no analytic for- 
mula for this extension. "A" shows the 1/kpa value at which the 
frequency shift changes sign, for C = 40. In the inset are plotted 
the differences between various curves and the bosonic mean-field 
value. 



so in both Fig. [T] and Fig. we plot C = 30 and C = 50 
for comparison. This value of C should be compared with an 
atomic BEC with large scattering length, for which C « 141 
& 

Monte-Carlo (MC) simulations l24l l25ll yield a curve for 
the radial breathing-mode frequencies significantly higher 
than the published experimental data fl2il . However, there are 
currently "word of mouth" indications that a new generation 
of breathing mode experiments from the Innsbruck group is 
in good agreement with MC data. If we use the MC results to 
deduce C, we obtain 90 < C < 120 (where our result is most 
highly constrained by the MC point [25] at 1/kpa = 1). By 
contrast, the range of acceptable values from the cloud size 
experiment is 40 < C < 65 (where our result is most highly 
constrained by the data point at B — 735G in Fig. 0. That 
there is no overlap might be of some concern. However, it 
can be shown that the MC data itself 1 25] is not entirely con- 
sistent (within error bars) with the cloud size data [8]. This 
is not unexpected since these Monte Carlo simulations, which 
are based on the fixed-node approximation |25] establish an 
upper bound to the equation of state; at the same time, it is 
quite plausible that these first generation cloud size experi- 
ments may be revised in due time. 

It should be stressed, finally, that the range of validity of our 
low-density formula is, strictly speaking, outside the range of 
fields accessed in the data which are closer to unitarity. Fu- 
ture experiments which penetrate more deeply into the low 
density regime will be needed to determine the parameter C 
in more detail. As a parenthetical observation, one can see by 
comparing both figures that the curves associated with differ- 
ent values of C in the cloud size experiments show a smaller 
deviation from each other, than do the curves associated with 
the collective modes. This is because the cloud size is propor- 
tional to the 1/5-th power of the dimer-dimer scattering length 
j3l and is thus less sensitive to the details of the equation of 
state. 

In conclusion the central contribution of our paper is the 
proposal for a new many-body ground state wave function of 
the atomic Fermi gas with BEC-BCS crossover. This wave 
function yields a ground state energy consistent with the cor- 
rect dimer-dimer scattering length add in the BEC regime |7]. 
It also reproduces the Lee- Yang term, whose numerical value 
is also fixed by add- We have also examined the next term 
beyond the Lee-Yang correction in a phenomenological fash- 
ion, building on cloud size data and collective mode experi- 
ments. One has to view this phenomenological analysis with 
some caution since experiments are performed close to uni- 
tarity whereas the theory presented here is most appropriate 
in the deep BEC. More experiments on the cloud size and 
collective mode frequencies in the deep BEC regime will be 
required to narrow down a range of the parameter C which ap- 
pears in our theory, and thus determine the appropriate equa- 
tion of state. In the present paper we have had to make due 
with experiments close to unitarity and with word of mouth 
information which indicates that the measured radial breath- 
ing mode frequencies are in a state of flux, and tending to now 



5 



yield good agreement with the Monte Carlo data. With that 
proviso the value of C can be as small as 40 and as large as 
120. Finally, we stress that there is a distinction here between 
a BEC of fermionic dimers and an atomic Bose gas with large 
scattering length where it has been shown 1 2 1 ] that in the latter 
Cps 141. 

Note Added- After this work was completed we became 
aware of a later related study by Levinsen and Gurarie 12611 
which produced the correct mean field term in the equation of 
state, but which did not address the Lee- Yang term, and the 
other corrections considered here. 

This work is partially supported by NSF under NSF- 
MRSEC Grant No. DMR-02 13745. We thank J. E. Thomas 
and C. Chin, as well as Q. Chen, E. Braaten, A. Bulgac, 
G. Shlyapnikov and D. S. Petrov for helpful communications. 
We thank anonymous referees of this manuscript for some 
suggestions. 
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